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We discuss the equilibrium condition for a liquid that partially wets a solid on the level of in- 
termolecular forces. Using a mean field continuum description, we generalize the capillary pressure 
from variation of the free energy and show at what length scale the equilibrium contact angle is 
selected. After recovering Young's law for homogeneous substrates, it is shown how hysteresis of the 
contact angle can be incorporated in a self-consistent fashion. In all cases the liquid-vapor interface 
takes a nontrivial shape, which is compared to models using a disjoining pressure. 



I. INTRODUCTION 

The equilibrium condition of a liquid that partially 
wets an homogeneous solid substrate has been addressed 
since Young [T], who found that there is a well-defined 
contact angle 9y- This angle minimizes the macroscopic 
free energy E of the liquid and provides the boundary 
condition for the free surface: 
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where 7, 7^; and 7^^ represent the surface tensions of the 
liquid- vapor, solid- liquid and solid- vapor interfaces re- 
spectively. The surface tensions are defined as the excess 
free energy per unit area, in the particular geometry of 
a planar interface between two unbounded phases [U [3] . 
When approaching the contact line where the three in- 
terfaces meet, the geometry changes dramatically and 
the force balance can no longer be expressed in terms of 
surface tensions. Instead, the interface deforms at small 
scales to establish a nontrivial equilibrium shape [H [5] . 
Macrocopically, however, one recovers a wedge of angle 

The small scale structure of the contact line is much 
less understood, but very relevant for problems such as 
line tension in nanofiuidics [TJ HI IHl HO] j or moving con- 
tact lines that are out-of-equilibrium down to molecular 
scales HH Ha 113 m [IS IS]- to avoid fitting parame- 
ters, these problems require an explicit treatment of the 
force balance within the range of molecular interactions. 
Another basic phenomenon not described by Young's law 
is contact angle hysteresis: chemical inhomogeneities or 
roughness of the solid substrate can trap the contact line 
in a potential well, allowing for a range of possible macro- 
scopic contact angles [H |T71 [TB] . A common interpreta- 
tion is to assume that the contact angle is selected at a 
scale smaller than the inhomogeneities, and to consider 
Young's law as a local boundary condition [121 1201 
Strictly speaking, however, this interpretation is not self- 
consistent because surface tensions are no longer well- 
defined at molecular scales. Alternatively, one could ar- 
gue that the global free energy E is dominated by con- 
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FIG. 1: Experimental verification of Young's law by |6]: the 
measured contact angle O^xp versus the prediction Oy from 
independent measurements of the surface tensions. 



tributions from the bulk of the drop, and described by 
an average over the inhomogeneities, jsv — Isi |22[ 123] . 
While this approach has been very successful in recent 
years for wetting on textured substrates [23], it predicts 
a unique value for the contact angle and hence does not 
capture the hysteresis. 

These observations raise the question of how the con- 
tact angle is selected, in particular at what length scale. 
Rather than imposing the angle as an external boundary 
condition at a large distance from the contact line [25], 
one would like to derive it directly from the force balance 
at molecular level. 

There are various approaches to address this problem 
within a continuum theory, but no general consensus has 
yet been achieved. A popular strategy consists of adding 
the microscopic effect of liquid film confinement (disjoin- 
ing pressure) and the macroscopic effect of interface cur- 
vature (Laplace pressure) [7] [TH [53] . A consequence of 
this approximation is that Young's angle is only recov- 
ered if a precursor film is present on the substrate. Con- 
ceptually, both the disjoining and Laplace pressures have 
the same physical origin: the force on the interface arises 
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because the geometry of the liquid deviates from a semi- 
infinite flat domain. A priori, these pressures need not be 
additive and a more rigorous approach consists of mini- 
mizing the total free energy of the system with respect 
to the interface shape h{x). From this, Getta and Di- 
etrich [8] derived a nonlocal integral equation for h{x) 
within Density Functional Theory in the presence of a 
mesoscopic precursor film. They demonstrated that the 
addition of Laplace and disjoining pressures is in fact a 
local approximation and they addressed the quantitative 
limitations of such a local approach. 

The situation without precursor film, typical of more 
hydrophobic substrates [27l [28], has been analyzed by 
a number of papers [5) UHl 133 131] • Surprisingly, these 
arrived at contradicting conclusions. Merchant and 
Keller [5] derived a nonlocal integral equation for the 
interface profile that has essentially the same structure 
as in [S] . From asymptotic analysis they showed that the 
profile macroscopically approaches Young's angle, while 
the inner structure of the contact line displayed nontrivial 
oscillations. On the other hand. Hocking [30j and more 
recently Wu and Wong [3F have identified an exact so- 
lution of a perfectly straight wedge down to molecular 
scale. Although unnoticed by these authors, the corre- 
sponding wedge angle is not Young's angle as obtained by 
Merchant and Keller, which raises an intriguing paradox. 

In this paper we revisit the equilibrium condition for 
the liquid- vapor interface at small length scales. We aim 
(i) to explicitly identify the connection between micro- 
scopic and macroscopic length scales, and thus the se- 
lection of the contact angle, (ii) to apply this to contact 
angle hysteresis, (iii) to clarify some of the mentioned 
contradictions in the literature. In view of the latter, 
the descriptions in the paper will in places be somewhat 
extensive. 

The analysis is restricted to mechanical equilibrium 
and treats the interface as infinitely thin, hence ignoring 
thermal fiuctuations and the corresponding details of the 
density profile. These simplifications allow writing down 
a free energy functional E[h] that contains all length 
scales but yet remains accessible for analysis. The capil- 
lary pressure 11 is generalized as the functional derivative 
of the free energy with respect to the interface shape h. 
n turns out to be exactly equal to the potential energy 
associated to the intermolecular forces. The equilibrium 
condition is then simply a constant capillary pressure, i.e. 
an iso-potential free surface: 11 = A. 

A first result is that the shape of the interface is never a 
perfectly straight wedge, due to regularization of the van 
der Waals interactions at a distance £, typically several 
Angstroms This regularization is crucial, as other- 
wise the total energy diverges and surface tensions are 
not defined. For homogeneous substrates it is found that 
the solution asymptotically approaches Young's contact 
angle in the large scale limit, effectively reproducing the 
result by Merchant and Keller [5] . The small scale struc- 
ture, however, does not exhibit the oscillations seen in [S]. 
We then identify how the large scale angle is selected in 



the neighbourhood of the contact line, at the scale £ of 
the regularization, and provide a self-consistent descrip- 
tion for contact angle hysteresis. Finally, we compare 
the generalized capillary pressure to the usual disjoining 
pressure. 



II. MACROSCOPIC ANALYSIS 

The benchmark for our analysis is the usual macro- 
scopic theory that we wish to recover from the inter- 
molecular forces. Here we present a formal derivation 
(similar to }33]), which can later on be generalized to in- 
clude microscopic interactions. We consider profiles that 
are translationally invariant in one direction, so that the 
problem reduces to finding h(x). The free energy of the 
liquid film reads (per unit length y), 

pb i*a poo 
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where from now on we write jiv = 7. The factor Vl + h''^ 
is required to compute the surface area (arclength) of the 
liquid- vapor interface. To minimize this energy under the 
constraint of a fixed liquid volume one needs to introduce 
a Lagrange multiplier A, and then consider the variation 
of £ = E ~ \ J dx h. 



5E = E[h + 5h]-E[h] 
= T{h)5h-T{a)5a + 
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Besides the usual partial integrations, we find contribu- 
tions da and Sb that correspond to contact line variations. 
These arise because there is no external force to constrain 
the contact line positions, so one can not impose 5h — 
at the boundary. 

At equilibrium this variation S£ — for arbitrary 
Sh(x). In the bulk of the liquid, the vanishing of the 
integral yields the Euler-Lagrange equation: 
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which for T(h, h') of ^ reduces to the usual Laplace 
pressure condition 

H = -7K = A, (6) 
where a — h" /{l+h'^)"^/^ is the curvature of the interface. 



3 



Similarly to the Hamiltonian in classical mechanics, 
this differential equation has a first integral, 
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that is conserved along the free surface. This can be seen 
by multiplying ^ by h' and integrating along x. To close 
the problem, we need to determine the value of tJmacro 
from the contributions at the contact line x = a, 6 in Q. 
The variation of the contact line da is not independent of 
6h(a). Through Taylor expansion one finds the relation 
Sa = —6h{a)/h'{a). As the contributions should vanish 
for arbitrary Sh, this provides so-called natural boundary 
conditions: 



at X = a.x 



(8) 
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Comparing to ([T]), one recognizes that these boundary 
conditions simply yield ^macro = 0. For the problem at 
hand we thus find 
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At the contact line this indeed gives Young's law Q. 



III. GENERALIZED EQUILIBRIUM 
CONDITIONS 

A. Formulation 

Before analyzing the microscopic free energy, let us 
reiterate the structure of the analysis. The bulk equilib- 
rium condition is obtained from the functional derivative, 
which basically corresponds to a constant pressure. We 
therefore introduce the generalized capillary pressure 11 
as the functional derivative, yielding the equilibrium con- 
dition in the bulk: 
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where A is, again, the Lagrange multiplier associated to 
incompressibility. 

The generalized invariant Q can be constructed from 
integration 



dxh'{x) [Tl{x) - A] 
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which from ( 10 1 vanishes for all values of x. In the macro- 
scopic calculation we closed the problem from the varia- 
tion of the contact line, which gave a term T{a)5a. This 
contribution is due to the discontinuous jump from r^ct 
to Fdiy across the contact line. We anticipate that on a 
molecular level such a discontinuity does not occur, as 
the relevant energies vary smoothly with the thickness h. 
This implies that in the microscopic model we lose the 
boundary condition. Instead, the solution will be selected 
internally from the balance of intermolecular forces. 



B. Interpretation 

The capillary pressure can be interpreted as a purely 
geometric effect. For a perfectly flat interface between 
two semi-infinite phases it is zero, because a virtual dis- 
placement of the interface will still yield a surface be- 
tween two unbounded phases. However, any deviation 
from this geometry will result into a force on the inter- 
face. The best known example is of course the Laplace 
pressure for curved interfaces, for which a virtual dis- 
placement leads to a change in surface area. Another 
example is the disjoining pressure for a thin horizontal 
liquid film on a solid substrate, whose thickness falls 
within the range of intermolecular forces. In this case 
the molecules near the surface feel a change of environ- 
ment, due to a replacement of liquid for solid molecules. 
We emphasise that both pressures have the same phys- 
ical origin and are captured within the single definition 
(flOl. 



The integral Q represents the horizontal component of 
the total force acting on the liquid between the contact 
line and the location x. Namely, the horizontal force per 
unit length y acting on a vertical slice of liquid reads: 



(ff = -U'{x)hdxdy 



(12) 



The horizontal component of the force is then obtained 
from integration by parts 

f- = - r diW{x)h{i) 

dy Ja 

= h{a)Il{a) - h{x)Il{x) + / dS; h' {S;)U{S:) 



- g. (13) 

As we have seen from the macroscopic analysis, this hor- 
izontal force balance should yield Young's law for ho- 
mogeneous substrates. In the presence of contact angle 
hysteresis, however, the invariant Q provides a nontrivial 
communication across length scales. We show below how 
it translates the microscopic force balance into a macro- 
scopic contact angle. 



C. Microscopic free energy and capillary pressure 

We consider pairwise molecular interactions, 0Q^(|r' — 
r|), where a and (3 can represent molecules in the liquid 
(1) or solid phase (s). We then follow [29j, by considering 
density-density correlations 



Pa/3(j',r') ~PaP/3 5a/3(|r'-r|), 



(14) 



where g^fs is the pair correlation function and we take pi^s 
constant over the liquid domain (C) and solid domain (5) 
respectively. This is equivalent to the so-called sharp- 
kink approximation in Density Functional Theory [8J . In 
the absence of external forces, this then yields the energy 
functional 
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E[h] = 



dv j^dv'MW-v\) 
dv f dr>„(|r'-r|) 



dr / dr'(bsi{\r' -r\), (15) 
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energy per unit volume at the free surface, due to the 
presence of liquid and solid molecules. The equilibrium 
condition is thus that h{x) is an equipotential [Sj [3D]. 
Mechanically, a gradient in potential energy would lead 
to fluid motion. Note that the recent paper [31^ uses the 
potential on the free surface to estimate the disjoining 
pressure, but then resides to a local approximation for 
the functional derivative, similar to equation (^. 



where we renamed 0Q/3(r) — gai3{f)4'ai3{T)- As this repre- 
sents the excess free energy, we subtracted the bulk liquid 
energy emerging from an infinite domain of interaction. 
The factors 1/2 arise since all pair wise integrations are 
counted twice by the double integrations. By definition, 
the gaseous domain is the complementary of the solid 
and liquid domains. One thus obtains an equivalent sys- 
tem if the gas is replaced by vacuum and if the gas-liquid 
interaction potential is subtracted from the liquid-liquid 
and the solid-liquid potentials. So, once the potentials 
are expressed in terms of surface tensions, one rigorously 
finds the same result if the gas phase is virtually replaced 
by a vacuum. 

We again consider profiles that are invariant along the 
y direction, bounded by contact lines at x = a and x — h. 
The domains of integration then become 
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The interface profile h{x) now appears explicitly as the 
boundary of the liquid domain. This makes it easy to 
evaluate the functional derivative, i.e. the capillary pres- 
sure, as 

n(a;) = / dx' / dz' (flu {x' — X, z' — h{x)) 
Ja Jq 

fOO j'OO 

dx' / dz' ifii {x' , z') 

-oo Jo 

/oo /"O 
dx' / dz' (fsi {x' — X, z' — h{x)) . 
-oo J —oo 

(17) 

Here we conveniently integrated out the invariant y di- 
rection, so that 



D. Surface tensions 

Once the molecular interactions have been specified, 
one can compute the surface tensions [H [3]. This is 
important in order to establish a connection with the 
macroscopic limit. An elegant way to obtain the liquid- 
vapor tension 7 is to directly derive the Laplace pres- 
sure from (17 1 in the macroscopic limit. In this case 



the boundary z ~ can effectively be replaced by —00, 
while a,b — ±00. A Taylor expansion h{x') — h(x) — 
h'(x)u + h" (x)u'^ /2 and changing variables u — x' ~ x 
then yields ^ IMj 
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with 



7 = -- y drr 4>ii{r). 



(19) 



(20) 



Indeed, this is precisely twice the energy required to sep- 
arate two semi-infinite liquid domains from contact to 
infinity. 

The other surface tensions can be derived from the 
disjoining pressure of a perfectly flat horizontal fllm, 
n(jisj(/i), through the connection pSl [31] 



1 + lsl- Isv 



dh ndisj(/i). 



(21) 



This relation expresses that ndisj(^) is the derivative of 
a surface free energy that has the correct macroscopic 
limits Fwct and Tdry. Taking h{x) = h{x') = h and 
a,b = ±00, we find upon integration 
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(22) 



Note that (20 1 is a special case of this result, where the 
solid is replaced by liquid so that ^sv = and jsv = 7- 
For a more rigorous thermodynamic treatment of surface 
tension we refer to [3|. 



(p{x,z) 



dy (f> [ \/ x'^ + y"^ + z^ 



(18) 



The pressure at the interface can thus be split into a 
part due to liquid-liquid interaction and a part due to 
solid-liquid interaction, i.e. 11 = 11;; -|- 11^;. Inspection of 
the integrals reveals that these are simply the potential 



IV. HOMOGENEOUS SUBSTRATES 
A. Asymptotic analysis: Young's law 

We now compute the equilibrium contact angle from 
the microscopic interactions, in the case where the solid 
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surface is perfectly homogeneous. In the macroscopic 
analysis of Sec. Young's law arises from variation of 
the contact line position, da. If we consider such a varia- 



portional to /J**""^ dz. Since h{a) = by construction, 
this variation does not give a nontrivial boundary condi- 
tion. As the contact angle should still emerge from the 
horizontal force balance, we evaluate the invariant Q by 



tion for the free energy (15 1, we find a contribution pro- 



integration (111. If the solid substrates are spatially ho- 



mogeneous, the solid- liquid potential is a function of the 
film thickness only: 



Tlsiih) 



dx' 



dz' ipsi {x', z'-h). (23) 



It contributes to the invariant as 

Qsiix)^ f ' \hU,iCh)- (24) 
Jo 

Integrating to a large distance from the contact line, h ^ 
£, this can be further reduced to 



Gsiioo) = 71" / drr (t>si{r) 
Jo 



(25) 



where we used the definition of surface tensions (|22| of 
the previous paragraph. 

The contribution of the liquid-liquid interaction is 
more subtle, as it involves the (unknown) equilibrium 
shape h{x), 



Guioo) 



h{x') 



dx h'{x) 
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dx' I dz' ifiii {x' — X, z' — h{x)) 

dx' I dz' ifu{x',z') , (26) 

where we took a = and b = oo. However, it was 
shown by Merchant and Keller [5 that the integral can 
be evaluated analytically if one assumes that for large h 
the interface approaches a wedge of a well-defined slope, 
h{x) ~ xtanOoo- The trick is to consider the difference 
between the real h{x) and the straight wedge. Merchant 
and Keller then showed that this difference vanishes due 
to the double integration over x and x', irrespective of 
the shape h{x). In Appendix [A| we clarify the crucial 
steps in their analysis. The end result depends only on 
the value of the asymptotic slope: 

CJ/;(oo) = 7 + 7cos6'oo. (27) 

Taking A = 0, the invariance G ^ Gu + Gsi — yields 

7 cos 6*00 + 7si - Isv = 0. (28) 

Comparing to ([T]), we indeed find 6rx! — Oy from Young's 
argument. So whatever the shape h(x), the equilibrium 



solution approaches Young's angle for h ^ £. For finite 
drop volumes of typical size R, corrections are of order 
Xh cx 'yh/R, so that the asymptotic analysis is justified 
for ^ < /i < i?. 

To recapitulate, the variation of the contact line po- 
sition docs not provide a boundary condition in the mi- 
croscopic theory. Instead, there is an internal selection 
of the solution for which one can derive the asymptotic 
angle 6^0 — Oy ■ This is a completely microscopic demon- 
stration of Young's law. Imposing a wedge angle differ- 
ent from 9y, hence G ^ 0, will force the solution to pass 
through a minimum or a singularity, as only the solution 
with G ^ can reach h = 0. 



B. Equilibrium profiles 

To analyze the inner structure of the contact line, we 
numerically solve the integral equation 11 = A, with the 
pressure taken from (17 1. We consider the following ef- 
fective interaction 



-Cap 



for r > i 
for r < £. 



(29) 



This corresponds to attractive Van der Waals interac- 
tions, regularized below r — i, combined with a flat pair 
correlation function, g(r) = 1. For simplicity we take 
the same regularization length £ for the liquid-liquid and 
solid-liquid interactions. Using Sec. HID we can readily 
compute Young's angle as [27] : 



1 - cos 6Iy = 2 1 



Cll 



(30) 



For details on the numerical algorithm we refer to Ap- 
pendix [c| 

Figure |2] shows an equilibrium drop profile, with 9y = 
0.7. Close to the contact line we observe a change of 
contact angle from a microscopic value, 9^, to the equi- 
librium angle, 9ac = 9y. Indeed, the dotted line indi- 



cates Young's angle computed from (30 1, confirming the 
asymptotic analysis. This equilibrium angle is attained 
at a typical height i (set to unity in all figures) above the 
solid. At larger distances, outside the range of the inter- 
actions, the profile is simply the cylindrical cap expected 
from macroscopic theory. 

In Fig. [3] we present the microscopic angle as obtained 
for various wetting conditions. The microscopic and 
macroscopic angles coincide only for 9y = 7r/2, for which 
the solution is a perfectly straight wedge. For more hy- 
drophilic drops, however, the microscopic angle is always 
much smaller than Young's angle. In Appendix [B] we 
derive an approximate relation based on local theories: 



1 



V (2 - cos 9y 



cos^ 9\ 



(31) 



see also the Discussion section. This provides a good 
qualitative description of the numerical results (dotted 
line, Fig.|3|. 
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FIG. 2; Equilibrium profile of a drop with interaction ( |29[ ), 
with 9y ~ 0.7. Lengths are expressed in terms of I. Top: at 
large scales, the shape is a cylindrical cap. Bottom: a zoom at 
the contact line region reveals a variation of the slope, from 
a micrscopic angle 6^ to the a macroscopic angle Ooo- The 
dotted line confirms that 9aa = Oy- 




FIG. 3: Relation between microscopic angle 
angle Oy 



and Young's 
Dotted line: 



Dots: numerical resolution for (291 
interpolation based on local approximations ( |31[ |. Solid line: 
~ dy, plotted for clarity only. Note that — 6y only at 
7r/2. 



Our numerical results gave no indication of the sur- 
prising oscillatory behavior of h{x), as observed by Mer- 
chant and Keller [S]. In their simulations, the macro- 
scopic wedge angle was imposed as a boundary condi- 
tion at infinity, which is different from our finite volume 
nanodroplets. Despite this technical difference, the inner 
structure of the contact line should be the same. Investi- 
gating the microscopic parameters in [5] it seems that the 
asymptotic boundary condition was not correctly taken 
at Young's angle [35]. This means that macroscopically 
t/ 7^ in their simulations. Enforcing a solution that 
at the same time has h = 0, where G = 0, the solution 



would have to pass through a singularity. Looking at 
the cusp-like structures in [5], we think this is a plausi- 
ble explanation. This problem could not possibly occur 
in our simulations, as the drop was free to establish the 
equilibrium shape. 



C. Paradox of the straight wedge 

The findings discussed above contradict the result by 
Hocking |30| . who analytically found equilibrium solu- 
tions that remain a perfectly straight wedge at all scales. 
Hocking's solutions exist for all wetting angles, while in 
our simulations the wedge was only found for Oy = 7r/2. 
The paper [30J starts with the same expression for the 



capillary pressure (17 1, and then uses unregularized Van 



der Waals interactions 



(t>al3{r) 



Cap 



for r > 0. 



Imposing the shape h(x) — xtanO this yields 



H. 



wedge 



6/^3 



Cll 



m 



(32) 



(33) 



with F{9) = 5 — I cos + \ cos^ 9. In order for the pres- 
sure to be constant along the interface, for all h, it is 
then argued that the equilibrium angle follows from 



Csl 
Cll 



1 



Cll 



1/4 



(34) 



where the latter approximation holds for small angles. 

Although the paper claims to be consistent with 
Young's law and this is not the case. Comparing 
(34 1 to (30 1, it is clear that 6h ^ Oy- Crucially, the 



interactions ( |32[ ) lack a small-scale regularization and di- 
verge when r — !■ 0. Integrating over this singularity gives 
a divergent contribution, so that the surface tensions are 
effectively infinite. As we have seen, the wedge solutions 
cease to exist when a regularization is introduced. 

The expression ( 33 ) nevertheless has a useful interpre- 



tation as the scaling for Van der Waals interactions when 
h ^ £. In fact, it has been used to derive the approxi- 
mation (31). For small contact angles, the contribution 
F{9) oc can be neglected so that the disjoining pres- 
sure is unaffected by the small inclination. For 9y = tt/2, 
however, the term in brackets vanishes and the effect of 
Van der Waals forces completely disappears. In this par- 
ticular case, the wedge indeed is an exact solution as 
observed in our simulations. 



V. CONTACT ANGLE HYSTERESIS 

A. Analysis 

The microscopic determination of 9oo follows from 
asymptotic analysis and does not depend on the homo- 
geneity of the solid surface. One can distinguish two 
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types of inhomogeneity that both lead to contact angle 
hysteresis: geometrical roughness of the substrate and 
variations in surface chemistry. In the latter case, the 
analysis of the Q integrals can be done explicitly. Let us 
therefore consider a spatially varying surface potential. 



-7(l + cos6ly)^'°(/i)-e7Cos(ga; + a)^'^(/i), (35) 



where we normalized dh'i>'^'^{h) = 1. Such a variation 
could reflect the crystalline structure of the solid, which 
gives rise to a spatial patterns of the solid density ps on 
molecular scale. The chosen form can be seen as a per- 
turbation expansion for small inhomogeneities or as part 
of a Fourier decomposition of the chemical variations. 

The liquid-liquid interaction is completely unaffected 
by this surface inhomogeneity, so we can evaluate the 
condition Q = 0, 



cos 9^ 



cos Hy 



(36) 



where 



dxh'{x)'^\h) cos{qx + a). (37) 



It is clear from (36 1 that the macroscopic angle 6*00 can 
take a range of values, and depends sensitively on the 
phase a at the contact line position (taken at a; = by 
definition) . 

The calculation of Gq{a) requires the slope h'{x) of 
the nontrivial equilibrium solution. One can, however, 
get an estimate for the hysteresis by introducing the ap- 
proximation h'(x) ~ tan^y: 

/•oo 

Gq{a) K / dh'i'^{h) cos{qh/ tandy + a) 
Jo 

= cos a 3? I ^ ' ^ 



where 



\tan ( 



*(g) 



sin a 3 ME' 



tan Oy 



dhe"'''^^{h), 



(38) 



(39) 



is the Fourier transform of ^'^(/i) times the Heaviside step 
function. The amplitude of the hysteresis, characterized 
by the advancing and receding contact angles 9a and Or, 
can then be estimated by 



cos 9a ~ cos 9y ~ e 
cos 9r ~ cos 9y + e 



tan^y 

q 

tan^y 



(40) 
(41) 



This result has two intuitive limits. Since the length 
scale appearing in ^'^(/i) will again be £, the hysteresis is 
controlled by the dimensionless wavelength 



Q 



q£ 



tan I 



(42) 



In the limit Q <C 1, one finds |^'| = 1, because of normal- 
ization. It this case the variation of the surface chemistry 
is slow and one can locally treat the surface as homoge- 
neous, with the relevant surface tensions varying between 
±e. This is the usual 'macroscopic' interpretation of hys- 
teresis based on surface tensions. In the opposite limit, 
Q 3> 1, one finds \^\ cx l/Q 0. This scaling is due 
to the singularity in the Heaviside function, which deter- 
mines the large scale behavior of the Fourier transform. 
In this limit, the spatial variation is so quick that it is 
averaged out and the hysteresis disappears altogether. 



B. An illustration 

The approximation above will now be compared to 
complete numerical solutions in the presence of hystere- 
sis (Fig. |4|i). We again consider the interaction (29 1, for 



which we can compute ^"(/i) for a homogeneous sub- 
strate. We consider spatial variations characterized by 
\l/i(/i) = ^'^{h), which gives (Appendix |c]) : 



2 

M 



{t/hf 
8~9{h/e) 



for h> I 
for;^<l (43) 



This function is indeed normalized to unity. The ampli- 
tude of the hysteresis can now be computed by evaluating 
the integral transform |^|, which is plotted in Fig. [ij; as 
a function of Q. 

We have numerically simulated droplets for this inter- 
action. We have imposed the drop center to coincide with 
a maximum the wall potential has a maximum, so that 
mirror symmetry yields two identical contact lines. The 
volume was kept constant, yielding a discrete set of solu- 
tions of varying contact angle 9oo ■ As can be seen from 
Fig. |4^, the contact line positions for these metastable 
solutions are separated by roughly one wavelength. The 
drops shown in this figure are the only ones we could 
obtained for this volume, reflecting the bounds of hys- 
teresis. 

As expected from our analysis, the precise phase a of 
the potential is slightly different for the various drops. 
This is shown in Fig. |4]3, in which the numerical results 
are compared to ( 38 1 . This approximation accounts quite 
well for the value of 9a, as well as for the phase shift 
with respect to T\si , whose extremal values are located at 
a = and 11. However, it signiflcantly underestimates 9^- 
We attribute this to the strong deviation from the wedge 
approximation near the contact line (c.f. the flattest drop 
in Fig. lit 



VI. DISCUSSION 

A. Internal selection of the contact angle 

We have shown how the macroscopic contact angle 
emerges from a microscopic force balance, within a mean 
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FIG. 4: (a) Metastable solutions for drops of equal volume in 
the presence of hysteresis. The phase of the interface poten- 
tial is sketched below (note that the solid-liquid interface is 
still taken perfectly flat). Simulation parameters: 9y ~ 7r/4, 
e = 0.427, 2 = 3. (b) Contact angle versus the phase of the 
potential Q. The numerical values (circles) and the approx- 
imation (38 1, (c) Absolute value of I*!'!, approximating the 



strenght of the hysteresis for wavelength Q. 



field continuum discription. There is a crucial role for 
the integral G, which is invariant along the entire in- 
terface. Physically, this invariant expresses horizontal 
force balance on an arbitrary cross-section of the liquid. 
The macroscopic angle can be computed by equating the 
value of Q for a macroscopic wedge to its value at the 
contact line. The height over which the integral attains 
its asymptotic value is simply £, the range of the interac- 
tions. In this respect, the contact angle is indeed deter- 
mined at a molecular scale. 

A striking feature of the microscopic model is that the 
entire solution follows from the pressure balance, without 
the necessity of an external boundary condition. Hence, 
there is an internal selection of the contact angle. This 
is very different from the usual macroscopic theory, for 
which the contact angle has to be imposed as a boundary 
condition. Besides this conceptual difference, the internal 
selection also works when surface tensions are no longer 
well-defined. This occurs for contact angle hysteresis in 
the case where variations of the surface chemistry are 
of the same scale as the range of interaction (expressed 
in Sec. |v] as Q ~ 1). Again using the invariant Q, the 
microscopic model naturally provides a range of contact 
angles. 

At this point we would like to stress that the free en- 
ergy analyzed in this paper oversimplifies the physics near 
a contact line: it ignores thermal fluctuations, while the 
continuum approach starts to break down at molecular 



scales . The present analysis should therefore be con- 
sidered as a model calculation, where one can explicitly 
identify the selection mechanism. It would be interesting 
to combine the present approach with molecular dynam- 
ics simulations, or more rigorous density functional calcu- 
lations to improve the physical reality. The results for the 
(submolecular) angle 6^, which resolved the paradox of 
the straight wedge solutions, should be interpreted within 
this context. Having said that, there exist experimental 
measurements of variations between macroscopic and mi- 
croscopic angles, for systems with sufficiently long ranged 
interactions [HI [in]- A very similar situation is encoun- 
tered in electrowetting, where an applied electric field 
induces a change of the macroscopic contact angle [37] . 
It was recently shown, however, that close to the contact 
line one recovers Young's angle [551155] . 



B. Local approximations 

We are now in the position to test local approxima- 
tions based on a free energy functional E = J dxT{h, h'). 
Such a formulation has the great practical advantage that 
the equilibrium condition reduces to a second order dif- 
ferential equation rather than a nonlocal integral equa- 
tion. The standard approach based on a disjoining pres- 
sure [3 [25] yields 6*^ = 0, independent of Oy- We refer to 
Appendix|B]for details. This is often interpreted in terms 
of a precursor film, which naturally has h' — 0, although 
strictly speaking this is not necessary. Comparing to nu- 
merical profiles of the nonlocal theory, we see that this 
provides a good approximation for small contact angles. 

The approximation clearly fails for large angles, for 
which the profile is much closer to a straight wedge, and 
0^ ^ 0. We therefore suggest an alternative local approx- 
imation: 

T{h, h') = 7 (v^l + h'^ - cos Oy) fih), (44) 

which has perfect wedge solutions, i.e. 9^ — Oy- We 
again refer to Appendix [B] for details. The numerical 
solutions obtained in Sec. IV B appear to be in between 
the straight wedge solutions and those from the standard 
models with disjoining pressure. The approximation for 
9^ in (31 1 has been obtained from interpolation between 



these two local models. The advantage of the form (44 1 



is that it provides an internal selection of the contact 
angle. This means that when applied to the nonequilib- 
rium situation of a moving contact line, the solution is 
self contained and no longer requires an additional con- 
dition borrowed from equilibrium. It remains to be in- 
vestigated how the local approximations compare to the 
nonlocal theory in the dynamical case. 

Let us conclude by noting that we have not been able 
to come up with a self-consistent derivation of a local 
approximation of the free energy for arbitrary 9. As sug- 
gested in [8], this is perhaps not possible due to the in- 
trinsically nonlocal character of the problem. 
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APPENDIX A: CALCULATION OF Gn 

We provide a graphical representation of the analy- 
sis by Merchant and Keller [5] to compute Gu- Expres- 
sion ( 26 1 can be seen as the integral over h, represented 



schematically by the black dot moving along the inter- 
face (Fig. ^p), of the potential energy at the interface 
evaluated at the position of the black dot. This potential 
itself is an integral over the entire liquid volume (light 
gray zone) of the intermolecular potential. Note that by 
considering ipu we already integrated out the invariant 
direction parallel to the contact line. Diagramatically, 
we have represented the potential as the interaction of 
the black dot with the dark grey rectangle, integrated 
over all positions of the open square. 

The trick is to compare the result for an arbitrary 
shape that tends to a given angle 9 (Fig. 5d), to that for 
a perfect liquid wedge of same angle (Fig. 5 The differ- 
nce between the two cases can be expressed as an integral 
over h (moving black dot) over the domain sketched in 
Fig. [5|i. This domain extends from the real interface 
(solid line) to the plane of angle 9 (dotted line). The dif- 
ference is counted positively if the interface is below the 
plane (light gray) and negatively otherwise (dark gray). 
If we interchange the position of the solid dot and the 
square dot, the potential energy due to the rectangle is 
exactly the same, but with opposite sign (Fig.|5^). Since 
both contributions are encountered through the double 
integration, the total difference vanishes by antisymme- 
try. This means that we can directly compute Gii by 
considering the simple wedge geometry, which gives the 
result ( [27] ) [5]. Note that this argument requires the as- 
sumption that the interaction depends only on the dis- 
tance between two points, but not on the spatial location 
of the points. If the interaction were inhomogeneous, the 
perfect antisymmetry would be broken. 



APPENDIX B: LOCAL APPROXIMATIONS 

The usual local approximation is to add the Laplace 
pressure and the disjoining pressure, giving the equilib- 
rium condition 





(c) 






1 






FIG. 5: (a) Qu can be interpreted as the total force on a liq- 
uid corner due to liquid-liquid interactions. We define h as 
the distance to the solid surface, (b) Graphical representa- 
tion of Qu as the potential at the interface at z = h, (black 
dot), integrated over all h (moving the black dots along the 
interface). The potential is decomposed as the contribution 
of the dark grey rectangular zone, summed over all positions 
of the square, (c) We compare the result for (b) to that for 
a liquid wedge of angle 6. (d) The differce beteen (b) and 
(c) is expressed by integration over a new domain (see text), 
consisting of positive (light gray) and negative (dark grey) 
contributions, (e) Inverting the two moving points, the con- 
tribution is identical but with opposite signs. 



This can be derived from the free energy 

Ti{h,h') = 7 (Vl + h'^ - l) +7(I-cos0y)/(/i), (B2) 

through the standard Euler-Lagrange description (see 
Sec. |ll]). The function f{h) is the integrated disjoining 
pressure, renormalized such that /(O) ~ 0, /(oo) — 1. 
The invariant corresponding to this free energy reads 



Gi^j 



( (l + fe,2)i/2 -l)+7(l-cosW(/^), (B3) 



so that the equilibrium condition Gi — yields 0p 



0, 



independent of 9y ■ 

For larger angles we suggest an alternative local ap- 
proximation: 



jK + Udisjih) = A. 



(Bl) 



T2iK h') = 7 ( Vl + - cos9y) fih), 



(B4) 
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which has a capillary pressure 



n2 - -iiifih) 



and a first integral 



dh V(l + /i'2)i/2 



cos Oy (B5) 



= 7 



(l + /l'2)l/2 



COS 6'y /(/l). (B6) 



The equilibrium solution has 6 — Oy for any value of 
/i, indeed representing a straight wedge. In particular 

The solution of the nonlocal equation appears to be 
bounded by the local approximations obtained from Fi 
and To estimate 0^, we therefore propose a linear 
interpolation 



= 1^0^,1 + (1 - Z/)0^,2 = {1-V) 



(B7) 



The weight v can be estimated from comparison to the 
pressure on a straight wedge of angle Oy, Eq. (33 1, as 
computed by Hocking [30) . Even though Hocking's wedge 
is not an equilibrium solution, it does provide the correct 
nonlocal asymptotic disjoining pressure for h ^ i. For 
the second local approximation one finds H2 = for a 
wedge oi9y, while the first approximation asymptotically 
gives (33 1 without the contribution F{9). We therefore 
equate 



1- 



Cll 



= ( 1 - ^ - F{e) 

Cll 



(B8) 



Combined with (B7l, this gives prediction (31) shown in 
Fig-i 



APPENDIX C: NUMERICAL SOLUTION OF 
DROP PROFILES 

This appendix provides details on the numerical resolu- 
tion of equilibrium drop shapes. We consider the effective 
interaction 



(i>al3{r) 



iov r > I 



(CI) 



-Cap/t for r<i. 
The constants Csi and cu can be related to the surface 



tensions, using (22 1: 



7(1 + cos 9y) 



37rC;; 

4^2 ■ 



(C2) 
(C3) 



1. Algorithm 



We perform the following iterative procedure for h{x) 
in order to find the equipotential surface. At step n, 



the (discretized) shape of the drop is noted hn{x). The 
potential corresponding to this shape, H„(a;), (see the 
second part of the appendix). The shape is evolved ac- 
cording to the equation: 



hn+i{x) = hn{x) + 5{<\ln> -H„(a;)) 



(C4) 



where < H„ > is the average potential. The parameter 5 
is fixed at a value sufficiently small to ensure numerical 
stability. The positions of the contact lines are evolved 
to ensure the mass conservation between steps: the drop 
spreads when the potential at the contact line is smaller 
than its average value; otherwise, it shrinks. 

In practice, we have observed that the result becomes 
independent of the discretisation when the mesh size be- 
comes lower than 0.05 {— 0.03 for the figures shown here). 
The convergence of the calculation, starting from a spher- 
ical cap takes approximately 300 steps (1000 for the fig- 
ures shown here). The spatial variations of the potential 
then reach the numerical noise level. 



2. Evaluation of the pressure 



We evaluate the pressure according to (17). The sec 



ond integral is constant and can thus be ignored. As the 
effective interaction is defined over two domains, we need 
to to separate zones of integration into those that have 
r > £ and r < £ respectively. First consider the solid- 
liquid interaction, i.e. the third integral of ( 17). If h > £, 
the integral is easily performed as 



lisiih) 



7r/2 p2n z'+oo 

d9 dip sin 9 (f) si ir)dr 

Jo J h/ cose 

m ■ 



(C5) 



For h < £, the domain has to be separated into various 
regions. Introducing cos a = h/£, we get for h < £: 



-2ttCs 



-27rCs/ / d9 sin 9 
Jo 



' h/ cos 9 



lisiih) - 

n/2 r-\-oc 

d9 sin 9 / dr/r'^ 

J h/ cos 6 

r^dr/£^+ [ dr/rA , 



which reduces to (for h < £), 



6£3 



(8 - 9ih/£) + 2{hl£f 



This potential H^i is indeed continuous at h 
obeys the property: 



(C6) 

(C7) 
and it 



dhUsiih) 



lsl-1- Isv = -7(1 + COS 6'y) 

(C8) 
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(a) 



J 




FIG. 6: (a) The interaction potential at the black point, due 
to the liquid zone at the same x as the open square (rectangu- 
lar zone), is decomposed into the two regions shown in light 
and dark gray, (b) The elementary configuration allowing to 
express the potential due to the liquid-liquid interaction is its 
contribution due to a band of liquid of height /i at a distance 
A from the point considered. 



Provided that we realise that "fn =0, the same formula 
can be apphed to the case of the liquid-hquid and sohd- 
solid interactions: 



cii 



(C9) 



The liquid-liquid potential can be expressed as: 



nb poo ph(x') 

Uu = dx' dy' / ~4>ii{r)dz' (CIO) 

Ja J-00 Jo 

with = {x' ~ 4- y'^ + {z' - h{x)Y . The integral 
over z can be decomposed into two contributions of the 
same form, as shown in Fig. [6]). Introducing the energy 
of interaction C(^i h) between a point and a rectangular 
zone of height h, infinite in the transverse direction, at a 
distance A, 

C(A, h) - p dy' p (VA2 + j/'2 + z'2) dz' 

(Cll) 

we get: 



Unix) = I dA (C(A, Hx)) + C(A, hix + A) - hix))) 

(C12) 



Importantly, C is defined as an odd function of h 

(C(A, — /i) = — C(A,— /i)) but an even function of A. 

In the situation where A > £, we introduce cos a = 
1/2 



cos^ 



1 



It simplifies into: 



and get: 



^<^^^U^^\..rcos^O^,^^^ (C13) 



df 



C ^ 2h^ + SA^ for A > ^ (C14) 

7 3A4(/i2 + A2)3/2 

Now for A < £, there is a circular screening zone of radius 
\/ P — A^. First situation, the screening zone is smaller 
than the height: 



167^2 



Stt ./o 



Tr/2 



d(p 



x6 A^da; ^1 1 A^dx 



For i"^ > A"^ > P — h? , it simplies into: 



(2/i2 + iA'^)eh 



2e 



4A2 



3A4(/i2 + A2)3/2 3A4 3£4 ^2 



X"^ 

(C15) 



(C16) 



Second situation, the screening zone radius \J P- ~ A2 is 
larger than the height h. We define cos (/Sq = hl\/P — A2 
and get: 



1 A^dx 




x^ 

'■^/'^ x6 AHx 1 A^dx^ 



A/a 



For t"^ — h? > A^, It finally simplies into: 



£6 



C _ 2h^P - A2 - /i2 / 1 ^ 4 (A2 + /l2) 
7 ~ 3(A2-t-/i2)7r l^A2 

2eh (3A2 + 2/i2) sin-i (VA2 + /i2) 



3A4 (A2 + /l2)3/2 



(C17) 



^ _ 3^ 2A^ 
37r V A4 ~ ^ ^ ^ 
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